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Feeding the black hole with condensing accretion flows: 
radiatively efficient and radiatively inefficient cases. 
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ABSTRACT 

We study the accretion flow of a hot gas captured by the black hole gravity in the 
presence of a thin cold accretion disk. Such geometrical arrangement is expected in 
Active Galactic Nuclei (AGN) and in galactic X-ray binary systems because both hot 
and cold gases are present in the black hole vicinity. Previous astrophysical literature 
concentrated on the evaporation of the cold disk in the classical heat conduction limit. 
Here we consider the inverse process, i.e. condensation of the hot gas onto the cold disk. 
We find two distinct condensation regimes, (i) In the classical thermal conduction limit, 
the radiative cooling in the hot gas itself force condensation above a certain critical 
accretion rate. Most of the flow energy in this case is re-emitted as X-ray radiation, 
(it) Below a certain minimum accretion rate, the hot electrons are collisionless and 
the classical heat flux description becomes invalid. We use the "non-local" heat flux 
approach borrowed from the terrestrial laser heated plasma experiments. Due to their 
very large mean free path, the hot particles penetrate deep into the cold disk where 
the radiative losses are signiflcant enough to enable condensation. In this case the hot 
flow energy is inconspicuously re-radiated by the transition layer in many UV and 
especially optical recombination lines (e.g., Lya, Ha, H/3) as well as via the optically 
thick disk emission. We derive an approximate analytical solution for the dynamics 
of the hot condensing flow. If the cold disk is inactive, i.e. accumulating mass for 
a future accretion outburst, then the two-phase flows appear radiatively inefficient. 
These condensing solutions may be relevant to Sgr A*, low luminosity AGN, and 
transient binary accreting systems in quiescence. 



1 INTRODUCTION 



An optically thick accretion disk (e.g., Shakura & Sunyaev 
1973) appears to be an important part of the accretion 
flow of gas into the black hole (BH) or a compact object. 
This conclusion is supported by spectral energy distribu- 
tions (SED; see Elvis et al. 1994; Ho 1999; and Fig. 1 in 
Gierlinski et al. 1999), double peaked emission line profiles 
(e.g. McClintock et al. 2003; and references in Ho 2003), and 
eclipse mapping of binary systems (e.g. Wood et al. 1986). 
Even for Sgr A*, a very dim source, there is now a suspicion 
that a cold inactive, i.e. not accreting, disk may be present 
(Nayakshin, Cuadra & Sunyaev 2004). 

There is also a significant amount of hot (T > 10^ K) 
gas at large distances from the BH in galactic nuclei. In the 
best studied cases this hot gas is observed as close as its cap- 
ture radius (e.g. in the giant eliptical galaxy and a LLAGN 
M87 [Di Matteo et al. 2003]; and in Sgr A* [Baganoff et al. 
2003]), meaning that accretion of this hot gas on the BH 
is unavoidable. In galactic binary systems, a diffuse hot gas 
may be present near the disk outer rim due to a shock in 
the hot spot, if the radiative cooling time is longer than 
dynamical time. A shocked captured stellar wind from the 
secondary is another source of hot gas (and it may also form 
the cold disk itself; Kolykhalov & Sunyaev 1980). 



Thus the accretion flow at large distances is often a two- 
flow problem (see Figure 0. Thermal conduction in a multi- 
phase medium leads to either evaporation of the cold gas, 
or vice versa, condensation of the hot gas (e.g., Zel'dovich 
& Pickel'ner 1969; Penston & Brown 1970). In reference to 
accretion flows, Meyer & Meyer-Hofmeister (1994; MMH94 
hereafter) were the first to study the mass exchange between 
the cold disk and the corona. Their pioneering study, ex- 
tended later by, e.g., Liu, Meyer & Meyer-Hofmeister (1997); 
DuUemond (1999); Rozahska & Czerny (2000a, b), showed 
that the cold disk evaporates at low coronal accretion rates. 
At large accretion rates the radiative cooling suppresses 
the corona (Rozahska & Czerny 2000a). The magnitude of 
the viscosity coefficient, a (Shakura & Sunyaev 1973), was 
shown to be extremely important. In all the models the evap- 
oration rate decreases very strongly with decreasing a (e.g. 
Meyer-Hofmeister & Meyer 2001). Recently, Liu, Meyer & 
Meyer-Hofmeister (2004) found condencing solutions for a 
sufficiently small value of the viscosity parameter. 

Here we aim to study the condensation process system- 
atically, i.e. attempting to cover a broad range of accretion 
rates in the hot flow. We first review the classical thermal 
conduction case. As expected, we recover qualitatively the 
results previously obtained by the above referenced authors. 
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Namely, when the heat flux is classical, we find condensa- 
tion at high and evaporation at low coronal accretion rates 
(32Jl. The critical accretion rate (at which no mass exchange 
between the disk and the hot corona takes place) is a strong 
function of a ( § & ^ . 

We then point out that for every value of a, there exists 
an accretion rate below which the classical thermal conduc- 
tion treatment becomes invalid since the hot electron mean 
free path, A, is long compared with the flow height scale, 
H. To study the problem in this limit, we use a modified 
(see §§ 12.11 fc 1^ prescription for the heat flux that is bor- 
rowed from the laser-heated plasma experiments. The heat 
flux is then proportional to the saturated heat flux coeffi- 
cient, (j) ^ 1. While this case still requires a future physical 
kinetics treatment, we show that there may exist an addi- 
tional condensation regime. In this regime the hot particles 
penetrate deep into the cold gas. The cold gas density turns 
out to be great enough to re-radiate the deposited energy 
away, enabling condensation. With the help of certain ap- 
proximations, we build an analytical solution (§§ for 
the radial structure of the hot flow in this case. 

We put our results into the context of the previous work 
in where we determine the type of solution for a given 
combination of parameters (accretion rate, radius, a and 
(j))- We discuss the expected spectra from the two types of 
the condensing flows in The classical thermal conduction 
condensing flow emits mainly in X-rays. Such solutions ap- 
ply at relatively high accretion rates, and therefore they are 
likely to be relevant to bright or medium-bright AGN, such 
as Seyfert Galaxies. Due to these high accretion/condesation 
rates, the cool disk is probably active, meaning that the SED 
of these sources would be dominated by the emission from 
the small scale accretion disk in the vicinity of the last stable 
orbit. Such flows would obviously be radiatively efficient. 

The second (non-local) condensation regime is applica- 
ble at lower accretion rates. It is therefore possible that ac- 
cretion disks fed by condensation at these low rates would be 
inactive, as in quiescent states of transient binary sources. In 
this case the hot gas, settling onto the cold disk, effectively 
stops accreting. The hot accretion flow is thus terminated at 
some large distance away from the central object, which im- 
plies that the radiation emitted by such a flow will be much 
less luminous than expected from a flow onto the black hole. 
Furthermore, these flows may be quite dim in X-rays since 
the raditive cooling in the corona itself is small. The cold 
disk serves as a huge cooling plate (radiator) for the hot 
flow. The SED of such flows would then be dominated by 
their infra-red/optical thermal-like bumps rather than UV 
or X-ray frequency regions. We suggest that this radiative 
"inefficiency" (rather a time-delayed accretion) of condens- 
ing flows may be one of the reason why some LLAGN, and 
the galactic black hole Sgr A*, appear underluminous. 



2 CLASSICAL THERMAL CONDUCTION IN 
TWO-PHASE FLOWS 

2.1 Basic physics of thermal conduction 

Figure Q shows a sketch of geometry that we envisage in 
this paper. The hot flow is sandwiching a very cold and thin 
disk located in the midplane. A thin (compared with R) 




R 



Figure 1. Geometry of the problem (shown for z > only). The 
hot flow (corona) sandwiches a cold accretion disk. A transition 
layer heated by hot coronal electrons (JSJ develops on the top 
of the cold disk. In the limit of coUisionless thermal conduction, 
both the transition layer and the disk are thick to coronal elec- 
trons, and are very geometrically thin. In contrast, the corona 
is thin to coronal electrons (since H ^ X), but is geometrically 
thick. Depending on conditions, the hot flow may be condensing 
(positive Vz) or evaporating. In the classical thermal conduction 
limit, the transition layer becomes optically thick to hot electrons. 
The layer's geometrical thickness also increases. 

transition layer develops between the corona^ and the disk. 
Since the cold disk viscous time is very long, the disk can 
be treated as stationary. 

The physics of the mass exchange between the two flows 
depends greatly on the heat flow from the hot to the cold 
gas. When collisions are abundant, i.e the mean free path for 
electron energy exchange. A, is smaller than the temperature 
scale height, Lt = T/|Vr|, the classical (Spitzer & Harm 
1953; Spitzer 1962) heat conductivity is used. The thermal 
heat flux, F, in this case is given by 

i? = Fci = -k{T)\7T , (1) 

where T is the temperature (we will assume that electron 
temperature is equal to the proton's temperature), k{T) is 
the conductivity coefficient given by 

fc(T) = 1.84 X 10~^r^''VlnAc erg cm^^ s"^ K^^ (2) 

where In Ac ~ 30 is the Coulomb logarithm. The heating or 
cooling rate due to thermal conduction is 

Ha^-divF. (3) 

For reference, 

A ~ lO^T^/n cm (4) 

(Spitzer 1962), where n is the gas density in cm"'^. If col- 
lisions are rare, i.e. A > Lt, the classical formula over- 
estimates the heat flux. The electron conduction is then ar- 
tificially fimited to (e.g. Parker 1963; Cowie & McKee 1977) 
the saturated heat finx: 

Fsat = 5(1) pel , (5) 

^ Throughout the paper, we shall use the terms "corona" and the 
"hot flow" interchangeably, although we do not assume that the 
origin of the corona is the cold disk. 



Condensing accretion flows 3 



where P and p are the gas pressure and density, respectively, 
= P/p is the local isothermal sound speed of the gas, and 
^ 1 is the saturation coefficient. Even though magnetic 
fields may be important in reducing the thermal conductiv- 
ity (see Appendix) , in what follows their effects are included 
only through the magnitude of <j). A more detailed discussion 
of thermal conductivity in the saturated regime is given in 



2.2 Classical evaporation or condensation 

Numerically accurate global evaporating solutions were ob- 
tained by MMH94, and that for condensation by Rozahska 
& Czerny (2000a). Here we only wish to establish the sign 
of the mass exchange between the two flows for given local 
conditions in the hot gas. For the sake of simplicity and gen- 
erality, we approximate the geometry (Fig.0 as planar, and 
neglect possible winds. In this setup the problem is similar 
to that of an infinite plane-parallel two-phase system stud- 
ied by Zel'dovich & Pikel'ner (1969), and Penston & Brown 
(1970). In such a system, for a subsonic flow, the entropy 
equation (for adiabatic index 7 = 5/3 and for Vz > for 
downward direction; Fig. is 

5fes dT dF 2 /fis 

- pvz — — + — =q+- HT)n , (6) 
Ip. az az 

where z is the vertical coordinate, q+ is the heating rate in 
erg/cm^/sec, A(T) is the usual optically thin cooling func- 
tion, and fcfl and p are the Bolzmann's constant and the 
mean molecular weight in the corona. Note that the subsonic 
flow assumption can be verified a posteriorly: as results of 
Cowie & McKee (1977) for evaporation of spherical clouds 
show, the flow is always subsonic when the heat flux is clas- 
sical. Also note that l?;^] <C Cs implies that the pressure is 
approximately constant throughout the thermal conduction 
interface. 

The boundary conditions for this equation are T{z — 
0) — Tc, F{0) = 0, and that the heating is equal to cooling: 
g+(0) - A(Tc)n2(0) = 0. Analogously, T(oo) = Th, etc. 

Multiplication of equation |S| by the heat flux and inte- 
gration from z = to z = 00 yields 

pv. = -| , (7) 
where 7i and I2 are the integrals 

/i= / 'drr^/^[g+ - A{T)n\z)] , (8) 



and 




respectively. Now, note that I2 is always positive, and hence 
the sign of the mass exchange depends only on that of the 
integral 7i. Further, the cooling rate is proportional to p^ 
while the heating rate often is proportional to just p, thus 
for coronal density ph high enough, the radiative cooling 
forces Ji < 0, and the hot flow will condense. Conversely, 
for low enough ph, the integral is positive and hence the 
evaporation is taking place. Evidently, for a given value of 
Th a unique value of pressure, at which the mass exchange 



is zero, may be defined. This pressure has been termed the 
"saturated vapor pressure" by Penston & Brown (1970). 

Consider now in more detail integral 7i . Above the peak 
at T ~ 10^ in the cooling curve, where A(r) ~ lO"^'^ erg 
cm"^ sec~^, following McKee & Cowie (1977), we approxi- 
mate the radiative cooling function as A(r) oc T'"'^ (until 
the free- free radiation losses become important). Since at a 
constant pressure n oc r"\ A{T)n^ oc T'^-^T''^ oc T"^ '', 
one could expect that the atomic cooling peak would sig- 
nfficantly contribute to the integral. However, the T^^^ 
factor strongly favors high temperatures, and as a re- 
suh we have / dTK{T){nTfT^''^ « dTA{T)T^^'^ oc 
JdTA(T)T^^'^ oc T°-^. This shows that the cooling contri- 
bution to Ii is strongly dominated by the highest tempera- 
tures, i.e. by T ^ Th- The heating rate is usually a weaker 
function of temperature and density than the cooling rate, 
thus the latter statement is even more appropriate. There- 
fore, approximately, we have 

h ~ T'J^iu - WWi^] « Q+~Q-, (10) 

where the vertical bar signifies averaging over the hot fiow, 
and the symbols Q+ and Q_ are the vertically integrated 
heating and cooling rates in the corona region. 

We shall recall that for typical heating and cooling 
mechanisms, a two-phase system without a thermal conduc- 
tion may be thermally unstable over a certain range in pres- 
sure (e.g. Field 1965; Krolik, McKee & Tarter 1983). In par- 
ticular, an initially uniform hot one phase medium will be- 
come thermally unstable and develop condensations when a 
critical pressure value is exceeded. Now, an important point 
is that this critical pressure and the saturated vapor pres- 
sure are usually rather close. For example, in the common 
Compton-bremsstrahlung case (McKee & Begelman 1990), 
the critical pressure is 0.95% of the pressure at which ther- 
mal instability in the hot gas starts to form clumps. The 
foregoing discussion and equation IIUI explain this result. As 
the integral 7i is dominated by temperatures T ~ T^, ther- 
mal conduction drives condensation at approximately same 
conditions as those when the thermal instability starts to 
form condensations in the hot fiow. 

In the context of an accretion flow, the critical pressure 
or density can be used to deflne a critical accretion rate, 
rhcrit. In this paper we measure accretion rates in units of 
the Eddington accretion rate, defined as LEdd/O.lc^, where 
7/Edd = 1-4 X 10'*^ Ms erg s~^, and Mg is the black hole mass, 
Mbh, in 10* Solar masses. In the Non-Radiative Accretion 
Flow theories (NRAFs, e.g. Narayan & Yi 1994; Blandford 
& Begelman 1999, BB99 hereafter), the radiative cooling 
term is not important {Q- ~ 0) as the gas is very tenu- 
ous. However, as the accretion rate increases, gas density 
increases, the radiative loss term eventually becomes im- 
portant, and a radiative collapse of the hot (one-phase) fiow 
occurs when Q- ~ Q+. The corresponding critical accretion 
rate, rhcrit ~ a^r^^'''^ (e.g. Narayan & Yi 1995; Esin 1997), 
where r4 is radius in units of 10'* i?s, with Rs = 2GMbh/c^, 
the Schwarzschild radius. The arguments made above for 
the classical thermal conduction show that the critical ac- 
cretion rate at which the mass exchange in the two-phase 
fiow is zero, rhdas, is quite close to rhcrit. In a simple model 
for the two-phase flow with 5+ — 9/2apCgf2K(l — T/Tvir), 
where the (1 — T/Tvir) approximately accounts for the effects 
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of winds, and free-free cooling, we found that the conduction 
drove evaporation at rhcias ~ 0.7mcrit- 

Summarizing these results, within the classical heat 
flux approximation, the hot flow evaporates the cold disk 
for accretion rates smaller than rhcias, and condenses for 
rh > rhcias- The classical condensation accretion rate rhdas 
is only slightly lower than the critical accretion rate merit at 
which the hot non-radiative flows without underlying disks 
would collapse to a thin disk. In the NRAF theories the ra- 
diative cooling term is not important, and hence these flows 
should drive an evaporation of the cold disk. 



3 CONDENSATION/EVAPORATION IN 
COLLISIONLESS LIMIT 

3.1 The non-local heat flux 

As accretion rate decreases, density decreases in non- 
radiative accretion flows. Thus at some point A (equation 
01 will become larger than the scale height H < R of the 
hot flow, and the classical heat flux description will become 
invalid (the exact condition for this is discussed in 
The common approach in the astrophysical literature is to 
limit the heat flux to the saturated heat flux value, as de- 
scribed in il2.1l However, although the saturated heat flux 
prescription remedies the situation as far as the maximum 
heat flow is concerned, it appears to be unreliable in re- 
gard to the correct estimate of the conductive heating rate. 
Namely, the latter is calculated as dFsa.t/dz = 5(jiPdCa/dz = 
54>P{cs/2T)dT/dz. While the heat flux is finite in the limit 
dT/dz — + oo, the thermal conduction heating is infinite. 
Physically, the maximum heating rate should be given by 
the heating rate of a cold test particle immersed in the hot 
gas (in fact it is far lower than that due to self-generated 
electric fields). 

A detailed calculation of the vertical structure of the 
transition region between the hot and the cold gas is be- 
yond the scope of the present paper (but will be taken up in 
Nayakshin & Sunyaev 2004, in preparation; NS04 hereafter) . 
The problem should be clearly treated with a physical ki- 
netic approach rather than a hydrodynamical one. However 
we can take guidance from the laser-heated plasma (LHP) 
research, in which the thermal conduction in the coUisionless 
(saturated) regime has been thoroughly studied via direct 
physical kinetics calculations and in-situ experimental mea- 
surements. At a physical level, the LHP research showed the 
importance of the high energy electrons that have very long 
mean free paths. These electrons carry most of the flux and 
their distribution function is not coupled to the the average 
quantities such as the mean temperature. 

Luciani, Mora & Virmont (1983; LMV hereafter), and 
Luciani, Mora & Pellat (1985) were able to formulate a "non- 
local" heat flux prescription that operates with hydrodynam- 
ical not physical kinetics quantities, and is yet able to cap- 
ture the majority of the relevant physics. The price that one 
has to pay for a greater self-consistency is that the heat flux 
becomes a non-local quantity. In particular, LMV expressed 
the heat flux, F^i, as an integral - a convolution of the clas- 
sical Spitzer's flux with a kernel that has a width of ^ A: 



where t(z, z') is the "electron optical depth" between loca- 
tion 2 and z' , 



X{z')niz') 



dz"n(z") 



(12) 



Fnl(2) 



dz' 



2X(z' 



- Fci(2:')exp [-t{z,z) 



(11) 



(Note that A used by LMV is about a factor of 2 larger 
than the electron equipartition A mentioned in H2.ll NS04). 
Formula llll has a simple physical interpretation: as the heat 
flux is dominated by electrons with energy E ^ bknT (e.g. 
see §IIA in Schurtz et al. 2000 for a nice account of this 
fact), their behavior is non-local. They are practically free 
streaming and thus the exponential factor in eouation llll is 
an attempt to model their transfer. 

Formula llll reduces to the classical limit when A <^ Lt, 
as expected. In addition, the non-local heat fiux may be 
shown to scale as the saturated heat flux'^ at an infinitely 
sharp temperature discontinuity, as it should. Most impor- 
tantly, the heating rate due to the non-local electron con- 
duction - the derivative of F^i (equation |^ - acts now on 
the kernel underneath the integral sign and hence it is finite. 

Although semi-empirical, formula |l II has been checked 
against numerical Fokker-Planck simulations (starting from 
Luciani et al. 1985). There have been improvements to the 
LMV kernel, and new non-local kernels were suggested (e.g. 
Albritton et al. 1986; also Epperlein & Short 1991, and refer- 
ences there). However they are all based on the same phys- 
ical ideas, and the LMV approach still remains the most 
popular choice (e.g. Ditmire et al. 1998; Schurtz, Nicolai & 
Busquet 2000) perhaps due to its relative simplicity. 



3.2 Vertical temperature profile 

The expected vertical temperature structure of the two- 
phase flow with a transition region is shown schematically 
m Figure H The a;-axis of the plot shows the column depth 
from the disk midplane, i.e. N = dz'n{z'). The corona is 
situated on the top and its temperature, Th ^ lO'^ K (may 
be much hotter, say 10® K). At the lower boundary of the 
hot flow, gas "temperature" rapidly decreases over a thin 
layer separating the corona from the "conduction heated 
zone" . In the both the classical and the saturated heat fiux 
theories, almost all of the heat carried by the hot electrons 
should be dissipated inside the region of the high temper- 
ature gradient. As the radiative losses there are small, this 
heating would be unbalanced and would lead to a vigorous 
evaporation. However, both the classical and the saturated 
heat fiux formulations fail to describe this region properly. 
The particle distribution function in the region is strongly 
non-Maxwellian as it consists of a mix of the hot and the 
cold particles. While the local "temperature" does define 
the mean particle kinetic energy, it is not related to the en- 
ergy dispersion at all. The non-local heat fiux predicts that 
only a small fraction of the heat flux, i.e. Az/A, where Az 
is the width of the transition zone, will be liberated in the 
region of the steep temperature gradient. Thus the evapo- 



^ The normalizing factor <p in equation |5] models complicated 
plasma collective effects, magnetic fields, etc., and may also be 



introduced in eauation llll 
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Schematic vertical structure of two-phase flow 
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Figure 2. Schematic vertical temperature profile of the two-phase 
flows (not to scale). The "cold disk" region is the bulk part of 
the disk. Directly above it is the layer heated and photo-ionized 
by emission of the gas above it (X-rays and possibly UV). The 
T ~ lO** — 10^ K region is conduction-heated by the energetic 
electrons penetrating directly from the hot flow. The depth of 
this layer is of order one mean free path of hot electrons. 



ration rate of this region is strongly reduced compared with 
the saturated heat flux predictions (NS04). 

The rest of the non-local heat flux enters the 
conduction-heated transition layer, in which the particle dis- 
tribution is dominated by the cold gas, and whose tempera- 
ture is 10^ K as this corresponds to the peak of the opti- 
cally thin cooling curve. Hydrogen in the layer is collisionally 
ionized. Below we shall show that the layer is dense enough 
to re-radiate away the deposited thermal conduction heat- 
ing. The emission occurs mostly in the optical and UV fre- 
quencies. The column depth of the conduction-heated zone 
is about one penetration depth, i.e. Np = An ~ lO^T^cm^. 
Since the hot flow by definition is optically thin to the hot 
electrons, the transition layer is actually optically thicker 
than the latter. 

Roughly a half of the radiation emitted by the 
conduction-heated zone, plus X-rays from the corona, is di- 
rected to and penetrate into a deeper cooler layer of the disk, 
photo-ionizing hydrogen there. This deeper layer is thus es- 
sentially an HII region. Finally, yet deeper is the main body 
of the disk. The disk is assumed to be optically thick and be 
at least as hot as the effective temperature corresponding to 
the heat flux from the hot flow {aBT^s — JLat). 



3.3 Energy balance in the conduction- heated zone 

The heating rate due to the non-local electrons penetrating 
the cold gas from the corona is 



dz 



= ^ exp l-riz)] = ^ exp hr(.)] , (13) 

^hc Ahc 



where we set the maximum of the non-local heat flux to 
equal to the saturated heat flux in the corona, at it should; 
t{z) is the optical depth of the hot electrons measured from 
the bottom of the corona, Ahc is the mean free path of hot 
electrons in cold gas, Ahc = Np/ric, and ric is the cold (that 
is local) gas number density. Note that the heating rate is 
simply given by the heat flux (JLat), divided by the effective 
distance to which the hot particles penetrate (Ahc), plus the 
attenuation factor (exp(— r(z)). 

Note that everywhere except the hot corona the pres- 
sure scale height is much larger than the temperature scale 
height, and hence the pressure is roughly constant with 
height (until the top coronal layers). With this, we have 
numerically 

-3/2 2 



(14) 



where Xe 1 is the fraction of free electrons in the transition 
layer, Ts = Tc/W^ K is its temperature, and Tg = Th/l(f 
K is the hot flow temperature. Introduction of x^^ above is 
a rudimentary attempt to account for the change in the en- 
ergy exchange rate, and is approximately valid as long as 
Xe Si [lnAc]~^. At lower Xe ionization of atomic Hydrogen, 
rather than Coulomb collisions, becomes the dominant heat- 
ing mechanism. 

The cooling rate of the transition layer is A{Tc)nc. The 
maximum of the cooling function for collisionally ionized 



gas is around Ai, 



10" 



erg/sec cm at temperature 



Tc ~ 10'' K (e.g. Raymond, Cox & Smith 1976). Comparing 
this with the heating rate given by equation 1141 we observe 
that this heating can be compensated by the atomic cooling 
ff Th ^ few X 10'' K or so. In addition, the non-equilibrium 
ionization state of the gas can further increase the cooling 
function by a factor of up to few tens (e.g., see Fig. 5 in 
Borkowski et al. 1990). This (i.e. A(T)n^ > Hni) is the 
regime that we are interested in. 

It is instructive to compare the cooling rates in the hot 
corona and in the conduction-heated layer. The former is 
~ 10~^^n| erg cm''/sec, whereas the latter is ~ lO"^'^^^ — 
■j^g-21 (j'^y'y^)2^2 ^ j g typically some 6-8 orders of magnitude 
larger. This is why the cold disk is such an efficient coolant 
for the hot coronal flow. 



3.4 Condensation of hot gas 

Up till now we have considered thermal conduction due to 
electrons only. The ion thermal conductivity is neglected 
in all three formulation - the classical, saturated and the 
non-local - for the reason that ions are much heavier and 
thus much slower than the electrons, thus carrying much 
less flux. However, in the saturated heat flux conditions, the 
electron heat flow is limited by factors of many to few tens 
compared with their free-streaming flux, most likely due to 
charge conservation and plasma instabilities arising when 
electrons stream past the ions. Thus the ion heat flow may 
become more important. 

Coulomb heating rates of cold gas by hot Meixwellian 
electrons and ions are well known (Spitzer 1962). In the 
problem of interest, the situation is far more complex since 
it is not well known whether and how the non-thermal ions 
will stream into the cold gas (e.g. see Balbus & McKee 1982). 
Below we shall make the simplest assumption, namely that 
the ratio of the non-local ion heating rate of the cold gas 
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to that done by the non-local electrons f equation I13II is the 
same as that for a Maxwellian distribution. Thus the total 
ion (assumed to be proton) and electron non-local heating 
rates are 



Htot = -ffni G{T,Th) 



(15) 



where the factor G is (Spitzer 1962; Balbus & McKee 1982) 
43 /p 



G = 1 + 



(16) 



(1-f 1836 T/rh)^/^ ' 

The factor /p < 1 is introduced to parameterize the fraction 
of the heat flux carried by the ions. Note that G — > 1 for 
T > Th/lOOO and that G^l+ 43 fp for T < Th/lOOO. As 
an example, for fp = 0.5, T = 10^ , and Th = 3 x 10^ K, 
G ~ 2.2, whereas for Tk = 10*, G = 5.5. 

Now, if the fraction of the non-local flux dissipated in 
the transition layer between the conduction-heated zone and 
the corona is negligible, and if h.{T)in? > //tot, condensation 
will prevail over evaporation. Since G is not much larger 
than unity for Th ^ 10* K, and since for larger Th the total 
heating rate decreases as T^^^^ (eauation I14II . it appears 
that the condensation will be taking place when the corona 
temperature is larger than the "condensation temperature" 

Tcond ~ 10^ K. 

It is likely that the hot flow temperature at the outer 
radius Ro (see H4.H is nearly virial. The latter exceeds Tcond 
for radii smaller than ~ lO^Rs, i.e. in most astrophysically 
interesting situations (e.g. cf. Fig. 5 in Narayan 2002). Also 
note that X-ray spectra of LLAGN are usually power-laws 
in 2-10 keV energy band, and thus on purely observational 
grounds the temperature of hot flows should be larger than 
10^ K for these. 



4 RADIAL DYNAMICAL EQUATIONS FOR 
HOT FLOW 

We now wish to write down simplified equations for the hot 
condensing flow with a goal of providing an example of the 
radial structure of such flows. We shall do this only for the 
newly found non-local condensation regime {^^, when the 
heat flux is saturated. 



4.1 Domain of applicability in radius 

We study the problem at radii smaller than R = Rc where 
Rc is the circularization radius of the hot gas. For one-phase 
accretion flows, the radial inflow velocity (e.g. see §5 Frank, 
King & Raine 2002) of the hot flow is vr ~ aCs{H/Rf ~ 
acs (for H/R ^ 1). Using this we obtain an estimate of the 
coronal gas density n at any radius R and a given accretion 
rate. We then compare the hot electron mean free path with 
R: 

4 -0.1 am" V-^/"* (17) 
H 

where m = M/M-Edd is the dimensionless accretion rate 
(Mfidd = iEdd/ec^, e = 0.1 is the efficiency of standard 
accretion, and L-^dd is the Eddington limit), is radius in 
units of 10*i?s, with Rs ~ 2GMbh/c^, the Schwarzschild ra- 
dius for BH of mass Mbh- Estimate 1171 was made as- 
suming that the gas is at virial temperature {Tvir{R) = 
GMBHmp/2kBR = 2.7 x 10^^"^ K; r = R/Rs). 



According to equation 1171 if A > i? at some radius R, 
then it is also true for all internal radii unless the accre- 
tion rate rh increases faster than oc r~^^'^ with decreasing 
R. Baring that, one can define saturation radius, -Rsat, via 
condition X/R — 1: 

2/3 



Rs^ 



2 X WRs 



(18) 



Below we shall explicitly assume that the heat flux is satu- 
rated. The domain of applicability of our solutions is thus 
R ^ Ro, where 



_Ro = min[_Rc, Rs^ 



(19) 



Formally, our treatment is valid only for non-relativistic 
electrons. However, for such high temperatures the electron 
mean free path in the cold gas is so large that its exact value 
does not particularly matter as long as most of the heat flux 
is liberated deep in the cold gas and is re-radiated away, 
allowing condensation of the hot gas. 



4.2 Geometry and assumptions 

The hot flow is sandwiching a very cold and thin disk lo- 
cated in the midplane (see Figure 0. As already discussed 
in 321 a very thin (compared with R) transition layer devel- 
ops between the corona and the conduction-heated layer of 
the disk. For our present purposes it is sufficient to approxi- 
mate both the disk and the layer as infinitely thin. Since the 
cold disk viscous time is very long, the disk can be treated 
as stationary. Under these assumptions the transition layer 
and the disk enter the dynamical equations for the hot flow 
as boundary conditions only. 

Due to thermal conduction into the disk, we expect the 
hot corona to be cooler than "normal" NRAFs that are ap- 
proximately at the virial temperatures far from the inner 
flow region (e.g. BB99; Narayan 2002). This makes the ther- 
mal pressure gradient not important in the radial momen- 
tum equation. The hot flow is thus in a roughly Keplerian 
rotation. Similarly, thermally driven winds are weak. 



4.3 Equations for the hot flow 

We now write height-integrated equations for a hot rotating 
flow with a possible mass exchange. Our equations are best 
understood in comparison with those of the well known so- 
lution obtained by Narayan & Yi (1994; hereafter NY94; see 
their eqs. 1-4), and the equations of MMH94. We begin with 
the mass conservation equation. In a one-phase accretion 
flow with no winds or mass deposition, the equation is writ- 
ten as d/dR [RHpvR] = 0, or M{R) = A-kRHpvr — const. 
Here M{R) is the radial accretion rate, vr is the radial flow 
velocity (vr > for accretion), p is the height-averaged den- 
sity in the hot flow. The vertical scale height H is introduced 
through the hydrostatic balance as H = Cs/Qk where Qk 
is the Keplerian angular velocity. Taking into account the 
exchange of mass in the vertical direction, the mass conser- 
vation equation is: 

^ [RHpvR] = pv^R , (20) 

where z is vertical coordinate with z — at the cold disk 
midplane. Since the radial pressure term is small, the radial 
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momentum equation (eq. 2 in NY94) is trivially satisfied. 
As the cold disk is also Keplerian, there is no exchange of 
specific angular momentum between the two flows and the 
angular momentum conservation equation (3 in NY94) is 
unaltered. With fl = we get 



M = AnRHpVR 



mK OR 



(21) 



The entropy equation should include (in addition to the 
usual terms) the thermal conduction flux, Fsat, and the hy- 
drodynamical flux of energy in the vertical direction. In the 



limit of subsonic vertical flows, c^, we write 



Q+ + pv^ 



5 2 GMbh 
2"' 2R 



(22) 



On the left hand side of the equation, Q+ is the vertically 
integrated viscous heating, and the second term is the me- 
chanical energy flow in the vertical direction. Q_ is the in- 
tegrated radiative cooling inside the corona: 



H 



(23) 



where ~ rrip, Th is the temperature of the hot flow, and 
A(r(i) is the usual optically thin cooling function. The cool- 
ing term due to the thermal conduction is just the saturated 
thermal conduction flux, -FLat- 

Our entropy equation 12211 is the height-integrated ver- 
sion of that in MMH94 (see their equation 8), with the fol- 
lowing exceptions. As explained in H4.2I we neglect winds 
here because we are mostly interested in cooler condensing 
solutions. Thus their side- way expansion term, important for 
z ^ R (see their eq. 5), is not included. Hence the last term 
on the right hand side of equation 8 in MMH94 is neglected. 
The radial entropy flow term (the third line in equation 8 in 
MMH94) is neglected for the same reason as winds. For Ke- 
plerian rotation, the conversion of gravitational energy into 
heat by a viscosity yields the following heating rate. 



Q+ = (9/2) apct 



(24) 



It is important to note that both the viscous heating and the 
thermal conduction cooling (equation]^ scale linearly with 
pCa- After some simple algebra and using GMbhH^ /2R? — 

ity is 



find from eg nation 1221 that the condensation veloc- 



M. 



MTV)pH 
+ 3^2 ci ^ 



(25) 



M 



(26) 



where A^c is the Mach number of the condensation velocity 
in the important limit of negligible radiative losses: 

3^ 2 

Note that had we not neglected radial advective cooling, 
then writing it as of Qadv ~ /advQ+, where /^dv ^ 1 is a 
parameter (NY94), the only necessary change to equation 
1^ would be to change a a{l — /adv)- 

We will be interested in the parameter space in which 
Mc > 0, i.e. condensation rather than evaporation. Never- 
theless, note that in the opposite limit eauations 1251 and 1261 
are not very accurate because we neglected winds. If a ^ </), 
the viscous heating dominates the energy balance, the flow 
is hotter, and it is unlikely that the winds can be rightfully 



neglected. The outflowing wind could simply carry the ex- 
tra heat away thus evaporating the hot flow itself (as in the 
BB99 solution) but not the cold disk. 



4.4 Boundary conditions 

At J? = Ro, we set gas temperature Th{Ro) to a given (^ 1) 
fraction of the virial temperature there. The density at -Ro is 
determined from a given M{Ro) value once vr{Ro) is found. 
Since Fsat oc p, the gas pressure, we set Fsat(2 — H) — in 
equation 1221 As we neglect the mass outflow from the top 
of the corona to infinity, we also write p{z = H) Vz{z = 
H) — 0. The value of the Vz{z = 0), i.e. at the bottom of the 
corona, is to be found from equations rather than prescribed. 



5 APPROXIMATE ANALYTICAL SOLUTION 

The first term in equation 1251 is a constant while the last 
one is proportional to density of the flow, becoming small 
at low accretion rates. We shall thus neglect this term in 
the non-radiative limit (see [jSJ. Inserting = A4cCs into 
equation 1201 and also substituting vr on its value found 
from equation I21II . we arrive at a second order differential 
equation that contains two unknown variables, p and c^: 



^ d 1 d 



fix 



PCs 



(27) 



This equation cannot be solved analytically or numer- 
ically in a general case. By introducing i;^ 7^ we added 
an extra variable to the accretion flow equations, and the 
number of independent equations is now smaller than the 
number of unknowns. This situation is well known in ana- 
lytical NRAF solutions where one has to introduce three free 
functions to obtain a solution (BB99). In our case there is 
no specific angular momentum exchange between the flows, 
and the heat flux is specified by -Fsat, therefore we need 
to parameterize "only" one physical quantity. We have cer- 
tain well motivated expectations about the radial temper- 
ature dependence. For example, for a pure non-radiative 
flow such as a NRAF with no thermal conduction cooling, 
Th{R) oc R~^ for large radii. On the other hand, thermal 
conduction tends to smooth out temperature gradients (for 
example within supernova remnants), and hence in the other 
extreme^, Th{R) ~ const. Thus, we set Cs = cq{Ro/RY, 
where ^ d ^5 1/2 and Co is the sound speed at Tio- 

If we define u = R/Ro, then equation 1271 can be 
re-written as 



pu 



3q -RoCo 9^ 



\pu 



AMc GMbh du'^ 
Finally, defining p = pu''^'^ , we obtain 



Here Ic is the dimensionless "condensation length" 



9a 



Cq-Rq 

AMc GMbh ^ 20(P - 18a GMbh 



3a Cq-Rq 



(28) 



(29) 



(30) 



3 In fact wc first studied the problem via numerical simulations 
in which we indeed observed this nearly constant Tf^{R) 
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We are mostly interested in the case (f) ^ a and there- 
fore we shall only explore the d — case below.'* In ad- 
dition, if a <C 0, the condensation length is small, i.e. 

< a/A4c ~ a/(j> <C 1. This circumstance facilitates find- 
ing an approximate analytical solution for equation 12911 : 



const exp 



1 

Lu 



(31) 

p/llu*' due to 



Indeed, d?p/du^ = p + l/llu^) 

the fact that 1/lcU > 1/lc ^ 1 for it < 1. 

Returning to the original variables, we explicitly write 
down the approximate solution: 



VR = 



3a eg 
21, Ro^k{Ro 

M{R) 



Mo — exp 




(32) 



Co — const <C Co , (33) 



(34) 



Here pc is the gas density at Ro, and Mo is the accretion rate 
at that point. Note that IcRq is roughly the radial distance 
over which the hot fiow will condense. At i?o — ^ ^ IcRo, 
there is little left of the hot flow. 

Note also that equation 13311 shows that the radial ve- 
locity is constant and is substantially smaller than the sound 
speed since we study cases when both a and TVfc ^ 1. Fur- 
ther, for Ic < 1/2 (recall that we assumed Ic <C 1), the 
accretion rate increases with R for all i? < _Ro, as it should 
for a condensing flow. Also note that the radial flow velocity 
is small compared with the sound speed and yet it is larger 
than expected for the usual one-phase hot flows for which 
VR,visc ~ aCs(H/R)^ < acs- 

Figure |3 shows M/Mq as given by equation 1341 for 
(j) = 0.2 and two different values of a, i.e. a — 0.02 and 
0.1 (solid and dotted curves, respectively). The accretion 
rates are normalized on their values at _R = _Ro- The hor- 
izontal dashed curve shows the case of no mass exchange. 
Note that for the case of very small a, m essentially plunges 
with radius, whereas for a = 0.1 (which is only a factor of 
2 smaller than </>), m ~ const for a range or radii before it 
reaches the exponential decline. 



5.1 Spectra from condensing flows 

For low luminosity systems, Comptonization in the corona 
is not important. Therefore we expect the coronal spectra to 
be dominated by the free-free emission as well as X-ray line 
emission (at the coronal ionization equilibrium) for Th ^ 10* 
K. The single temperature free-free spectrum is thus the 
hardest one that can be expected here; the photon spectral 
indices are then F ^ 1.4 or so. The spectra are softer for 
lower Th and also in the case when there is a broad distri- 
bution of temperatures in the corona. 

The maximum transition layer temperature is Tc ^ 10''' 
K. Gas at this temperature is a very efficient emitter of UV 



^ Nevertheless we note that (i) approximate solutions may be 
obtained for a general value of d, and (ii) for d = 1/2 there is a 
scale-free solution. 
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Figure 3. Radial accretion rate profiles in the hot flow normal- 
ized on Mo = M{Ro) (for _Ro = 3 X lO^Rs). The saturated 
heat flux coeflicient is <j> = 0.2 for all the cases. The M/Mg = 1 
(dashed) curve shows the case when the hot flow neither evapo- 
rates not condenses, which occurs when a ~ ij>. For a smaller than 
(f), the hot flow is condensing. When a is only slightly lower than 
<j), Ic ~ 1, and the condensation is slow (cf. the dotted curve). For 
a <C 0, condensation length Ic <^ 1, and the hot flow very quickly 
settles down on the cold disk. This latter case is illustrated by 
the solid curve for which a = 0.02. Finally, for a ^ <j>, the hot 
flow (but not necessarily the cold disk!) evaporates, outflowing to 
inflnity as in the BB99 solutions. 



lines, e.g. that of He^, Oxygen and other elements. The 
column depth and thus the emission measure of the gas at 
T ~ 10^ K depends on the exact value of the non-local ther- 
mal conduction heating, //tot (eauation llSH . which depends 
on the hot flow temperature mainly. If the temperature is 
below ~ 10* K, then the ~ 10^ K gas will re-radiate a sig- 
niflcant fraction of the condensing flow energy. In this case 
there should be a strong UV bump in the SED of the source, 
along with another bump of comparable magnitude at the 
effective temperature of the disk, which is in the infra-red 
to the optical frequency range depending on the parameters 
(black hole mass, condensation radius, condensation rate). 

If the hot flow temperature is greater than 10* K or so 
at the point where most of the condensation takes place, on 
the other hand, then the situation is qualitatively different. 
Since //tot is lower, the layer with T ~ 10^ K exists only due 
to local thermal conduction and is much thinner than in the 
opposite case (NS04). Bulk of the conduction-heated zone is 
at temperature T ^ 10'* K. Strong Hydrogen Lya, Ha, H/3, 
etc., line radiation is expected from these layers rather than 
UV emission. Thus in this case the UV region in the SED 
of the condensing flow is weak. 

As shown in Figure |5| below the conduction-heated 
layer there will be also a photo-ionized layer, an analog of an 
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Figure 4. The accretion rate, rhsat, below which the heat flux 
in the hot flow is saturated (coUisionless), as a function of the 
saturation parameter 0. AU the three curves are for a = 0.1 but 
varying radii as labelled above the respective curves. The rhsat 
curve is rather approximate in the region (p ~ a where the two 
different limits meet (see equation ■ Note that for the largest 
radius, the hot flows are saturated only for extremely small values 
of rh. On the contrary, for the smallest radius, the hot flows are in 
the coUisionless regime even for relatively large m unless (/> ^ 1. 



HII region. This deeper layer is expected to be a strong emit- 
ter of optical Unes, especially Ha, H/3. Therefore note that 
in the present model the vertical structure of the transition 
layer contains both coUisionally and photo-ionized emission 
regions whose relative contribution to the optical spectrum 
will vary depending on conditions. 

Finally, below the optically-emitting zone there will be 
the optically thick cold disk. Its effective temperature will 
be controlled by both the external energy deposition (the 
condensing flow plus starlight, etc.) and the internal energy 
dissipation. For LLAGN, this typically yields Tea ~ 10^ K. 



6 PARAMETER SPACE 
6.1 Saturated accretion rate 

First we specify the conditions under which the heat carry- 
ing electrons in the hot flow are coUisionless, i.e. the heat 
flux is saturated. The estimate is done at radius i?o, where 
we assume that the hot flow temperature is virial (in Hfci.4l we 
discuss R < Ro). Note that for a.(j) < 1 one has to reduce the 
electron mean free path A (equation |1J to A' = ^A. Indeed, 
the saturation coefficient (p is an ad hoc attempt (see Ap- 
pendix A) to account for the reduction of the electron heat 
flux due to magnetic fields, plasma instabilities, etc., and to 
be consistent one should also reduce the electron effective 
mean free path. The critical gas density can then be related 
to the accretion rate given that we know the radial velocity 
of the flow, vr. For (j) ^ a we use equation whereas for 
the case of <C a we can simply neglect the thermal con- 
ductivity altogether and use the usual viscous radial velocity 



scaling (see the first paragraph in H4.1II . Thus we arrive at 
the following condition for the thermal conductivity of the 
hot flow to be coUisionless: 



m < TTlsat = 



0.1 x/oATT <;/> r4 ii (j> > a 

0.1 a (j}r~^^'^, ii(j><a. 



(35) 



Figure^Jshows rhsat for three different radii and a same value 
for a- viscosity parameter, a = 0.1. The sharp drop in the 
curves around (/> ~ q is the artifact of our assumption tj)^ a 
that we used to build an approximate analytical solution in 
^ in the reality a numerical integration of the equations 
would yield a smoother transition between the two cases 
given in equation 1351 For radii as large as lO^iis, the hot 
Uow must be accreting at a dimensionless rate smaller than 
~ W'^ff for the heat Uux to be saturated. However the Uow 
is almost always saturated for "small" radii such as 10^-Rs. 



6.2 Radiative condensation 

As noted in spontaneous condensation may be tak- 

ing place even in an originally one-phase system without 
any thermal conduction if the in-situ radiative cooling is 
too large. Hot accretion Uows suffer a catastrophic collapse 
when the cooling exceeds the viscous heating. For a NRAF, 
this occurs at accretion rates exceeding rficrit ~ a^r^^^^ (e.g. 
Narayan & Yi 1995; although exact results depend on the 
cooling processes and geometry of the flow: Esin 1997). In 
the presence of the cold accretion disk below the hot flow, 
the additional thermal conductive cooling causes the latter 
to condense at a smaller pressure. In H2.2I we showed that 
this shifts the value of the critical radiative condensation 
rate rhrad downwards by a factor of order unity only, i.e. 
rii-rad ~ nT-crit (but See H6.4I belowl . 

Now, if m < msat, i.e. the heat Uux is coUisionless, and 
the thermal conductivity is large (<^ ^ a), we use the solu- 
tions obtained in JS] The radiative cooling contribution to 
the condensation can be neglected when K{T)pH/ (3/i^cf ) ^ 
Mc (see equation I26II . Evaluating this relation at R = Ro, 
assuming the free-free cooUng, A(r) = 2.4 x lO'^'^T^^^, one 
Unds the accretion rate, mrad, at which radiative losses start 
to dominate the energy balance in the hot condensing flow: 



nirad 



c 1/2 , ,3/2 -1/2 



(36) 



If m < rhsat but <C a, then the thermal conduction can 
be neglected and we again estimate that rhrad ~ ci^'r^^^^, 
as in a pure NRAF case. Combining together these results, 
we obtain the accretion rate above which condensation is 
driven by the in-situ radiative losses in the hot Uow: 



■Tlrad — 



2„-l/2 



if, 

u, 



and Fci > Fsat 
or Fci < Fsat 



(37) 



6.3 Division of the parameter space 

We can now divide the parameter space into different regions 
by the type of the Uow that occurs there. The division is 
governed by the following principles: 

• m < rhsat, rh < rhrad, ^ a: Non-Radiative Condensa- 
tion (NRC for short), the solutions obtained in this paper. 
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Figure 5. Left panel: The types of solutions occurring in the two-phase accretion flows at radius R = 10'' Rs and for a "large" value of 
a-parameter [a = 0.2), as a function of the dimensionless accretion rate and the saturation parameter <f>. The signs + or — reflects the 
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• rh < rhsat, rh < rhrad, ^ a: Non-Radiative Accre- 
tion Flow (NRAF), with a weak (since is small) collision- 
less condensation. The flow structure is similar to the usual 
NRAFs except the flow should be somewhat cooler. 

• m > rhsat, rn < mrad Classical evaporation of the disk 
by the hot flow (e.g. MMH94). "Classical" reflects the fact 
that the heat flux is in the classical, collision-dominated, 
regime. 

• m > rhrad: Independently of other parameters, the hot 
corona is condensing onto the cold disk because of a too 
strong in-situ radiative cooling. 

To illustrate the resulting division in the m— (j> parame- 
ter space, flrst consider large radii, such as i? = 10^ Rs, and 
large a (Figure]^ Left panel). For each of the segments in 
the parameter space, the respective type of solution is in- 
dicated, and the signs or "— " indicate the sign of the 
mass exchange between the cold disk and the hot flow; 
stands for the evaporation of the cold disk. 

At the largest accretion rates (above rh ~ 0.01), the flow 
condenses due to the radiative losses in the hot flow itself. 
Slightly below that there is a narrow zone of the classical 
condensation, which is enabled by the radiative losses inside 
the hot flow and the transition region. 

At yet lower accretion rates, the radiative cooling in 
the hot flow and the transition layer becomes too weak. The 
excess energy is used to evaporate the cold disk, as first 
discussed by MMH94. The hot flow also drives a thermal 
wind to infinity (MMH94), and in many respects is similar 
to the NRAF solutions (e.g. BB99; Narayan 2002), although 
the two-temperature assumption typical for NRAF is not 
required for these evaporating flows. 

Next, below m = rrisat, lie the parameter space where 
the heat flux is collisionless. As we found in this paper, as 
long as Th Tcond , the heating of the transition layer by the 
hot electrons can be counter-balanced by the local radiative 
losses, and thus the hot flow is condensing. When <j) ^ a, 



the approximate solutions developed in JS] apply. The hot 
flow is condensing due to radiative cooling in the transition 
layer. Unlike the classical condensation regime, however, the 
transition layer is very much cooler in this case (T ^ 10^ K; 
see I^J, the hot electrons are in a free-streaming colli- 
sionless regime. The hot flow can be completely quenched 
(terminated) in this regime as M(J? — > 0) ^ 0. 

Finally, at ^ a and m < rhsat , we have an analog of 
a hot non-radiative corona that we simply call a NRAF (al- 
though once again, it does not have to be a two-temperature 
flow as such). While the hot flow is still condensing in this 
case, the process is slow since <j) < a, and thus the hot flow 
will not be terminated as in the NRC case. Therefore, the 
NRAF (-) case may be quite analogous to the usual NRAFs 
except for the additional (weak) cooling and condensation 
provided by the thermal conduction. 

Consider now smaller radii ( e.g., R — 10'' i?s in Figure 
inj. The hot flow is now much hotter, and thus the column 
depth penetrated by hot electrons is much larger. There- 
fore the flow is much more likely to be in the collisionless 
regime. For sufficiently small radii, the hot flow can be both 
collisionless and radiative along the line separating the ra- 
diative condensation and the NRC solutions in the upper 
right corner of both panels in Figure |S| This means that for 
such small radii the classical evaporation regime completely 
disappears for large values of 4>. As can be seen in Figure |H| 
for (f) ^ a, the hot flow condenses for any rh, either due to 
in-situ radiative losses (high rh) or due to radiative losses in 
the transition layer (low m) as discussed in §§ 

6.4 On global nature of realistic two-phase flows 

As stated in the beginning of ^ so far we have only consid- 
ered the parameter space for condensation or evaporation at 
the outer radius of the hot flow, R = Ro, where we assumed 
Th — Tvir. However at smaller radii, where the radial (in- 
flow) velocity significantly decreases compared to the free- 
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fall value, one expects that the thermal conduction will be 
much more efficient in cooling the hot flow to Tu <C TVir- 

Let us determine how the value of rfirad changes in this 
case. The integrated heating rate, Q+, is the rate at which 
the inflowing gas liberates its gravitational energy, and hence 
it only depends on the accretion rate, m, and is indepen- 
dent of the coronal temperature. The radiative cooling rate, 
however, scales as Q_ ~ K(Th)p^ H . Since 
for accretion flows, and vr ~ aQ.KR{H / Rf' , we have p oc 
m(T^ir/Thf^^. Taking A{Th) ~ const, we see that the ra- 
diative cooling rate approximately equals the heating rate 
at 



rp . 5/2 

J-h \ 



-1/2 ( Th \ 



(38) 



which may be significantly smaller than that given by equa- 
tion |37| Evidently, the infiuence of the thermal conduction 
on the collapse of a hot fiow is very significant in this ex- 
ample in the cumulative sense. While a hot one-phase ac- 
cretion flow with Th — Tvir may be thermally stable to col- 
lapse (spontaneous condensations), the two-phase hot flow 
with same accretion rate would be much cooler, <^ Tvir, 
much denser, and would thus collapse onto the cold disk. 
The detailed calculations by Liu et al. (2004) show that this 
is exactly the case. These authors flnd a condensing solu- 
tion for a = 0.1 and accretion rate m ~ 3 x 10~^, which 
would be somewhat too low for the classical condensation. 
However their vertical temperature profiles show that the 
coronal temperature is only Th ~ 0.3rvir, making rh < riirad 
according to equation 1381 and enabling condensation. 

Even more complex situation may exist when the gas 
is evaporating at large radii but condenses at small radii, 
or the flow is in the classical thermal conduction limit at 
large radii and is yet in the saturated conduction limit at 
small -R. This complexity of the two-phase flows with ther- 
mal conduction is not something new, however, as accretion 
flow equations are well known to allow more than one self- 



consistent solution at same M and radius (e.g. Abramowicz 
et al. 1995). 



7 APPEARANCE OF CONDENSING FLOWS 

From the disk ionization instability theory (for reviews see 
Cannizzo 1998, Lasota 2001), it is well known that "cold" 
accretion disks may be in two states: (i) active, that is ac- 
creting, in which case the disk is approximately in a steady 
state, with the accretion rate equal to a constant indepen- 
dent of radius, (ii) quiescent or inactive, when most of the 
hydrogen is neutral. The disk then accumulates the mass 
and the accretion rate strongly decreases with decreasing 
radius R. 

In the case of an active disk, we expect that a steady- 
state will be reached in which the mass deposited by the con- 
densing flow is accreted with time onto the black hole. The 
luminosity of the inner accretion flow is then much larger 
than that of the condensing flow at large radii. Accordingly, 
the latter is a minor detail in the overall spectrum of the 
source. This situation is more likely to realize for the clas- 
sical thermal condensation that takes place at a relatively 
high accretion rate than the non-radiative condensation. 

In the opposite case of a quiescent disk, the situation is 
completely different since then the luminosity of the system 
is much smaller and the condensing flow may be the pri- 
mary radiation source. This astrophysically interesting case 
is more likely for the non-radiative condensation, and is dis- 
cussed below. 



7.1 Underluminous accretion flows 

If the cold disk is inactive, i.e. barely accreting, in analogy 
to quiescent disks in binary transient systems (e.g. see Can- 
nizzo 1998; Lasota 2001), then the mass condensing on it 
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may be stored there for very long times. Indeed, the viscous 
time scale at radius R is 



tvis 



3 X 10 years 



(39) 



where Ms is the black hole mass in units of 10* years, 0:2 
and T3 are the cold disk viscosity and the temperature in 
units of 0.01 and 10^ Kelvin, respectively. For example, for 
Sgr A*, with Ms ~ 0.03 and with Ts ~ 0.1 (e.g. Nayakshin 
et al. 2004), this corresponds to 10* years, far longer than 
the age of the bright young stars estimated at ^ 10^ years 
(e.g. Ghez et al. 2003). For a black hole mass of ~ 10 M©, 
tvisc ~ 30 years, a time relatively long for binary systems. 

For inactive disks, the accretion rate through the inner 
boundary. Mi , is much smaller than that at large radii (e.g. 
Nayakshin & Svensson 2001), Mq. Some rough estimates are 
possible. Suppose the condensing gas is deposited onto a ring 
with radius R = Ro, and that the initial mass of the cold 
disk, Md{t ~ 0), is negligible small. Thus the disk mass at 
time t is Md{t) ~ tMo. The time needed for the matter to 
diffuse to the innermost stable orbit is exactly the viscous 
time. The rate at which the black hole consumes the disk 
can be estimated as (see, e.g., §5 of Frank et al. 2002)^ 



M,(t)~^exp 



t 



Mo 



t 

ivisi 



■ exp 



.(40) 



This equation shows that as long as t ^ tvisc, the innermost 
accretion rate is very small compared with that of the hot 
condensing flow. The expected luminosity of the innermost 
disk is Li ~ Q.lMiC? . Comparing this with the luminosity of 
the condensing flow at i? = Rq, Lcond ~ (3i?s/ft))0.1c^ Mq, 



U{t) ^ Mit) 



Ro 
3Rs 



Rq 
3R^ 



t 

tvis 



exp 



tv 



t 



(41) 



Even if Ro 3> 3Rs, the luminosity of the inner disk is small 
as long as t ^ tvisc. 

Therefore, we see that the hot condensing flow in coupe 
with an inactive cold disk are "underluminous" when com- 
pared with accretion flows that carry the gas directly into 
the black hole (i.e. assuming the standard 10% efficiency 
of the mass to radiation conversion for accretion). The ef- 
fect is of course a time dependent one; upon averaging over 
t 3> tvisc, the average luminosity of the two- phase flow equals 
that of the standard accretion flow (Shakura & Sunyaev 
1973). 



T.3 SED for non-radiative condensation 

In the two-phase flows studied here, most of the thermal en- 
ergy resides in the hot corona. However, in the non-radiative 
case, the corona is cooled by thermal conduction instead of 
radiation. Owing to this, the radiation output of the con- 
densing flow is dominated by the radiation of the under- 
lying disk. X-ray emission for such flows should always be 
week compared with the thermal-like bump in the infrared 
to UV frequency region (depending on the black hole mass 
and other parameters). 



7.4 SED for radiative condensation 

If the hot gas condensation is initiated by the radiative losses 
in the corona itself (©, then most of the hot flow energy 
is emitted in X-ray frequencies. Baring the case of a very 
optically thin cold accretion disk, about half of this radiation 
is absorbed in the disk and is re-emitted in the infra-red 
and/or optical frequencies. Thus in this case the SED of 
the source should contain a similar amount of power in the 
X-ray and the disk components. 

Clearly, a case intermediate between this and that con- 
sidered in ii7.3l exists. when the hot flow persists for a decade 
or more in radius before succumbing to the condensation 
process (see end of tl6.41 . 



7.5 Broad double-peaked emission lines 

As we discussed in > I5.1I the gas sandwiched between the 
corona and the cold disk is a strong emitter of optical and 
UV (if hot enough) lines. X-rays from the corona also con- 
tribute to radiative ionization and emission of optical/UV 
lines. Also note that the hot completely ionized coronal gas, 
settling down on the disk, will have to recombine, again pro- 
ducing emission lines. Therefore a generic prediction of our 
model is the presence of strong emission lines that should 
be shifted and broadened due to disk rotation, resulting in 
double-peaked line profiles. We expect the lines to trace a 
broad range of conditions and hence be a mix of photo- 
ionized and coUisionally ionized cases. 



7.2 Missing or truncated inner disks 

If the cold disk is truly inactive, only the area of the disk 
actually covered by the hot flow will be emitting. If one 
estimates the disk extent via the velocity width of Hq line 
or via the SED of the source, then it may appear that there 
is no disk at small radii because most of the hot flow does not 
reach there. While it is quite possible that the inner disk is 
missing altogether, it may also be in the inactive state where 
it emits very little to be detectable, as explained above. 



^ Exact results will depend on the run of the cold disk temper- 
ature with radius, but the dominant exponential factor will be 
always similar 



7.6 Fluorescent Fe Kq lines 

Here we have considered the hot X-ray emitting flow at large 
{R ~ 10"^ — lO^iis) distances from the black hole. If such a 
flow condenses and liberates most of its energy at these large 
distances, when clearly any Fe Ka line due to reprocessing 
of the radiation in the cold disk would be "narrow", in ac- 
cord with the observations (e.g. Ptak et al. 2004; Dewangan 
et al. 2004). Furthermore, because the inactive disks may be 
powered by condensation of the hot gas rather than by the 
internal disk accretion, the column density of these may be 
much lower than that deduced from the cold disk equations 
with the accretion rate tuned to match the optical/UV out- 
put of the source. Thus the inactive disks studied here may 
actually be Thomson-thin, producing no detectable Fe Ka 
line at all. 
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8 CONCLUSIONS 

Here we studied the physics by which a hot flow above a cold 
accretion disk could condense. Such a condensation process 
is of a large practical importance for accretion in AGN and 
X-ray binaries. Without means to cool, the hot gas is an ex- 
tremely inefficient kind of fuel for the black hole (e.g. BB99) . 
Not only the gas radiates little, it also accretes little (for re- 
cent numerical simulations see Proga & Begelman 2003). In 
contrast, if the hot gas condenses onto a cold accretion disk, 
the gas may temporarily loose its viscosity, i.e. the ability to 
accrete efficiently, but it remains tightly bound to the black 
hole. With time, when accretion through the cold disk is 
restarted, the hot gas will accrete onto the black hole. Thus 
our condensing scenario for the accretion of the hot gas onto 
the black hole is much more efficient in feeding the hot gas 
into the black hole than the non-radiative hot flows (BB99, 
Narayan 2002). 

We have found two distinct condensation regimes: (i) 
radiative or classical, which takes place when the radiative 
cooling term in the corona is comparable with the viscous 
heating term (see Sj^J ^-nd the thermal conductivity is clas- 
sical; and (ii) non-radiative. When the flow density is very 
low and the radiative losses are negligible. However the hot 
particles have very long mean free paths that allow them 
to penetrate the cold gas directly. Their energy is then re- 
radiated by the dense cold layers (see The cold disk 
thus serves as a cooling plate, radiator, for the hot flow. 
This second condensation regime is not expected based on 
the previous literature that employed the classical or the 
saturated heat flux formulations. 

We also presented a simple analytical solution describ- 
ing the radial structure of the hot flow for the non-radiative 
condensation. The solution is obtained under the assump- 
tion of hot flow temperature constant with radius and is 
meant to be an example only. The detailed structure of the 
flow is not as important for the flow energetics as long as 
one can clearly deflne a radius R' where most of the corona 
condenses. As the hot gas condenses, the gravitational en- 
ergy released per unit mass of hot gas is ~ GMbh/2R' . If 
the cold disk is inactive, then most of the condensed mass 
is stored in the disk until the disk becomes massive enough 
(e.g. Siemiginowska, Czerny, & Kostyunin 1996; but note 
that galaxy mergers, etc., may be another "direct" way of 
triggering accretion outbursts in the inner parsec from the 
BH). Thus the luminosity of the condensing flow (while the 
disk is inactive) is quite small compared with that expected 
for standard accretion flows f UT.lH . 

Nuclei of LLAGN is an example where the two-flow ge- 
ometry is natural - e.g. see Fig. 2 of Ho (2003) and note 
that in reality the hot gas probably fills the whole available 
space because it is too hot to be confined (BB99). Currently, 
the bump- like infra-red feature seen in the SED of LLAGN is 
thought to be due to the accretion through the disk. Assum- 
ing that the flow in the disk is continuous down to the last 
stable orbit however leads to a paradox since then one would 
expect a far brighter source. In contrast, if the cold disk is 
powered by the condensing corona and there is no accretion 
in the disk itself, there is then no problem with the source 
being too dim. The current "non-radiative" condensation 
mechanism is thus somewhat inconspicuous in that the ra- 
diative output of such a flow is dominated by the emission 



of a warm atomic and ionized gas and not by the energy 
source - the hot gas. 

A similar picture may be relevant for quiescent disks in 
transient binary systems if not all of the accretion stream 
radiatively condenses in the hot spot (the hot spot emission 
is quite weak in many systems; see references in Lasota 2001, 
his §7). Thus in some systems the "disk" continuum and line 
emission may be excited by shocks in the hot spot, whereas 
in others (with a weaker hot spot), the emission may be 
produced by the diffuse hot flow condensing onto the disk. 

Since the accretion flow may span a broad range of radii, 
it is possible that different types of mass exchange solutions 
realize at different radii. For example, at a high enough ac- 
cretion rate, the hot flow may be condensing radiatively at 
large radii, whereas at small radii, where the accretion rate is 
signiflcantly reduced, it may condense by the non-radiative 
mechanism. It is equally possible for the accretion flow to 
be evaporating at large radii and condensing at small ones. 
Therefore we should keep in mind that in general a mix 
of the solutions is possible. This "unfortunate" multi-phase 
complexity of accretion flows is not something unique to 
the accretion process; after all the interstellar medium has 
a multi-phase structure too. 

Finally note that due to processes not taken into ac- 
count in the standard disk instability model (e.g. Cannizzo 
1998), for example photo-ionization of the disk upper lay- 
ers by starlight and X-rays; stellar impacts (Nayakshin et 
al. 2004), etc., there will always be a finite (weak) accre- 
tion onto the BH. This weaker accretion may power the jet 
emission. 

Our results also indicate that a left alone cold disk will 
not necessarily have to evaporate into a hot corona at a low 
m via the classical thermal conductivity (MMH94). Indeed, 
when a hot corona just starts developing above the disk, the 
coronal accretion rate initially is very low and it is thus in 
the coUisional, not classical, regime. 



9 ACKNOWLEDGMENTS 

We thank Rashid Sunyaev for his careful reading of an ear- 
lier draft of the paper, and the many useful suggestions and 
comments that he made that substantially improved this 
work. We also acknowledge very useful discussions with F. 
Meyer, E. Meyer-Hofmeister, B. Liu, H. Spruit & K. DuUe- 
mond. Finally, C. McKee is thanked for several illuminating 
discussions during the author's visit to UC-Berkeley. 



APPENDIX A: NOTES ON COMPLICATIONS 
NOT CONSIDERED HERE 

Magnetic fields Magnetic fields are not explicitly intro- 
duced in this paper. The reason is pragmatic - magnetic 
fields in accretion flows are known to be turbulent and time- 
dependent (e.g. Balbus & Hawley 1991) and there seems 
to be no hope of obtaining an analytical solution such as 
that found in ^ when thermal conduction and fields are 
both considered. The presence of even a very weak (far sub- 
equipartition) magnetic field strongly reduces the thermal 
conductivity in direction perpendicular to the field lines. 
The usual approach is to set the conductivity to zero in that 
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direction (e.g., see Balbus 1986). Thus, the effective value of 
(j) will be reduced by any non-zero magnetic field, but an 
essential point is that this reduction is a function of the field 
topology but not the field strength. The main unknowns in 
this regard are the field correlation length and the angular 
distribution of the field directions near the interface between 
the hot and the cold flows. 

The value of the effective viscosity coefficient a is also 
a function of the magnetic field, but a different function. 
While the correlation length and the angular directions are 
again important, the absolute magnitude of a is nearly lin- 
early proportional to the magnetic field pressure (e.g. Balbus 
& Hawley 1991). Furthermore, the effective value of a should 
be obtained by averaging over the whole hot flow, whereas 
for (f> only the intermediate layers are important. 

For these reasons there is a broad spectrum of theoret- 
ical possibilities in the magnitude of the ratio (j)/a. For a 
small correlation length and near equipartition fields, one 
expects (f) <^ a, whereas for weaker fields and a larger cor- 
relation length the reverse should be true, iji 2> a. 

One-temperature assumption We assumed here that 
the electron and proton temperatures are equal. This is 
clearly so at the cold disk and the transition layer. For the 
hot corona, the Coulomb interactions may be inefficient in 
keeping the two temperatures equal. However, even under 
the assumption that almost all the heating goes into pro- 
tons, and only the Coulomb interactions operate, the two 
temperatures in a NRAF are approximately equal at radii 
R > 10^ Rs (e.g. see Fig. 5 in Narayan 2002). Thus in the 
most "interesting" for us range of radii the one-temperature 
assumption should be quite reasonable. 

In the innermost regions of the flow the protons may be 
much hotter than the electrons. Clearly the heat flux is then 
dominated by the protons. However, this heat flux should 
be of order 4'pc'^, i.e. only slightly lower. Hence our main 
argument remains the same: saturated thermal conduction 
will induce condensation of the hot flow when the viscous 
heating is weak. The physics of cooling in the transition 
layer will however be different for small radii (e.g. Spruit & 
Haardt 2000). 
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